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SPECIFICATION 
OPTICAL PULSE RECONSTRUCTION FROM SONOGRAM 

RELATED APPPLICATION 

This application claims the benefit of U.S. provisional application Serial No. 
60/272,888, filed March 2, 2001 , the benefit of which is hereby claimed under 35 U.S.C 
119. 

FIELD OF THE INVENTION 

This invention relates to an optical pulse reconstruction from sonogram, and more 
particularly, to a method for measuring the optical pulse from its sonogram and an 
optical sampling system employing the same. 

BACKGROUND OF THE INVENTION 

Frequency-resolved optical gating (FROG) is most commonly used to measure the 
amplitude and phase of ultra-short optical pulses. In the FROG system, we measure the 
spectrogram, which is the field spectrum of an optical pulse under test temporarily 
gated by itself. Nonlinear optical materials are employed for such optical gating. 
Though the window function for optical gating is unknown, the amplitude and phase of 
the pulse are retrieved from the measured spectrogram by an iterative minimization 
algorithm. 

An alternative approach is the sonogram characterization method. In this 
measurement, after a pulse is frequency-filtered, the intensity waveform of the filtered 
pulse is measured by a cross-correlator which is based on optical mixing using 
nonlinear optical materials or two-photon absorption in photodiodes and 
semiconductor lasers. It is shown in D. T. Reid ("Algorithm for complete and rapid 
retrieved of ultrashort pulse amplitude and phase from a sonogram," IEEE J. Quantum 
Electron, vol. 35, pp. 1584-1589, Nov. 1999) that an iterative algorithm similar to that 
used in the FROG system, in which the window function for frequency filtering is 
assumed to be unknown, can be employed for pulse reconstruction from the sonogram. 

According to the algorithm proposed in above-mentioned article, however, 
time-consuming iterative calculations are indispensable for pulse reconstruction from 
the sonogram. 
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An object of the present invention is to retrieve the amplitude and phase of an 
optical pulse from its sonogram without iterative calculations. 

Another object of the present invention is to provide rapid pulse retrieval from the 
sonogram. 

Still another object of the present invention is to enable us to discuss the sampling 
pulse width required to reconstruct the pulse accurately. 

Still another object of the present invention is to provide a formula for 
accomplishing the above-mentioned objects. 

SUMMARY OF THE INVENTION 

According to the present invention, there is provided a method for measuring an 
optical pulse which comprises: filtering an optical pulse to obtain a frequency-filtered 
pulse, a transferor window function for said frequency filtering being given; measuring 
a sonogram, which is defined as the intensity waveform of said frequency-filtered 
pulse, to obtain a measured sonogram; and reconstructing said optical pulse by using 
said measured sonogram and said transfer or window function. 

The present invention also provides a formula for retrieving the amplitude and 
phase of an optical pulse from its sonogram. When the transfer function of the 
frequency filter is known, the pulse amplitude and phase are completely retrieved from 
the sonogram without iterative calculations by derived formula. The pulse 
reconstruction formula is practically important for rapid pulse retrieval from the 
sonogram. More importantly, it enables us to discuss the sampling pulse width required 
to reconstruct the pulse accurately. 

The present invention also relates to an optical sampling system including the 
sonogram characterization function. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing aspects and many of the attendant advantages of this invention will 
become more readily appreciated as the same become better understood by reference to 
the following detailed descriptions, when taken in conjunction with the accompanying 
drawings, wherein: 

Figure 1 shows experimental setups for measuring the sonogram in which the 
sampling pulse, which is synchronized with the pulse under test and has a width 
narrower than that of the pulse under test, is prepared; 



3 



Figure 2 shows experimental setups for measuring the sonogram in which the 
sampling pulse is the same as the signal pulse under test; 

Figure 3 shows a modified process for pulse reconstruction from the sonogram 
obtained in Fig.2; 

Figure 4A shows optical sampling systems including the sonogram 
characterization function in which the sampling pulse is obtained by compressing the 
signal pulse under test; 

Figure 4B shows optical sampling systems including the sonogram 
characterization function in which the sampling pulse is incident on the device under 
test (DUT), and the impulse response of the DUT is characterized from the sonogram; 

Figure 5 shows experimental setups for the optical sampling system having the 
sonogram characterization function; 

Figure 6 shows the measured transfer function of the bandpass filter for frequency 
gating; 

Figure 7 shows the sonogram of the output pulse from the bandpass filter under 

test; 

Figure 8 shows the intensity waveform and phase of the output pulse reconstructed 
from the sonogram; 

Figure 9A shows the auto-correlation trace and spectrum in which circles show 
those calculated from the reconstructed pulse, and solid curves are directly measured 
ones; and 

Figure 9B shows the auto-correlation trace of the signal pulse in which circles 
show those calculated from the reconstructed pulse, and solid curves are directly 
measured ones. 

DETAILED DESCRD7TION OF THE PREFERRED EMBODIMENT 

I. SONOGRAM MEASUREMENT 

In the sonogram measurement, after a pulse is frequency-filtered, the intensity 
waveform of the filtered pulse is measured by a cross-correlator which is based on 
optical mixing using nonlinear optical materials or two-photon absorption in 
photodiodes and semiconductor lasers. 

In the sonogram measurement, we have freedom to choose a sampling pulse width 
for cross-correlation. Two extreme cases are shown in Figs. 1 and 2. In Fig. 1, we 
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prepare a sampling pulse, which has a pulse width narrower than that of the pulse under 
test and is highly synchronized with the pulse under test. After the pulse under test is 
frequency-filtered, it is cross-correlated with the sampling pulse in order to measure the 
intensity waveform of the frequency-filtered pulse. In this case, we can determine the 
sonogram precisely as long as the sampling pulse width is short enough. 

On the other hand, in Fig. 2, the pulse under test is divided into two replicas. One 
of the replicas is frequency-filtered and then cross-correlated with the other. However, 
the sonogram thus obtained is not an actual one because the temporal resolution is 
limited by the shape of the pulse under test. 

n. PULSE RECONSTRUCTION FORMULA 

We first derive a formula for retrieving the pulse under test from its sonogram. In 
contrast to the sonogram characterization of the prior art, we assume that the window 
function for frequency filtering is given. In such a case, the pulse amplitude and phase 
are completely retrieved from the sonogram without iterative calculations by using the 
derived formula. 

Let the complex amplitude of the signal pulse under test be s( x ). When the signal 
pulse is frequency-filtered by a band pass filter whose transfer function is H{ a> ), the 
complex amplitude of the output pulse is given as 

s^t) = — L= / exp(jio't)S{u/)H(ui -uj')dto', (1) 
x/2ir J 

where S( co ) is the Fourier transform of s( t ). The sonogram is defied as the intensity 
waveform of the frequency-filtered pulse: 

G(u;.t) = \s,M)\ 2 . (2) 

We next discuss how we retrieve s( x ) from G( co , t), closely following the 
method described in L.Cohen, "Time-Frequency Distributions- A Review," Proc. 
IEEE.,vo\. 11, no.7, pp941-981, 1989. 

The characteristic function M(6 , x ) of the sonogram G(co,t) is defied as 

M(0.t) = J J G(uj,t)exp(j8t + jroj)dtdu. (3) 
On the other hand, the ambiguity function for the signal is defined as 
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A s {0, r) = j s~(t - ^r)s{t + ^r) cx V (j0t)di. (A) 
If the inverse Fourier transform of H{ go ) is h{t), given as 

Hi) = ^J ^Hexp(jw<)rfw, (5) 
the ambiguity function of /z(t) is similarly expressed as 

A h {8.r)= j h*{t- l -T)h{t+ l -T)ex-p{j6t)dt. (6) 

Then, the characteristic function M(6 , r ), defined by (3) can be expressed in terms of 
these ambiguity functions as 

M{B,r) = A s {B,T)A h {-6.r). (7) 
From (4), we have 

s"(t - ^r)s(t + ir) = ~ j A a (O,r)exp(-j0t)dO. (8) 

By letting t= % /2 in (8) and substituting (7) into (8), we obtain the following pulse 
reconstruction formula: 

If MiO.t) 

We find that when the transfer function H( co ) of the filter is given, the complex 
amplitude s(t) of the pulse is completely retrieved from the measured sonogram G(o>,f) 
by using (3),(5),(6), and (9). 

On the other hand, Reid discusses an algorithm for pulse reconstruction from the 
sonogram based on iterative calculations, where it is assumed that s(t) and H(co) are 
unknown. The pulse reconstructed from the sonogram by this algorithm should be the 
same as that given by (9). However, once the transfer function of the filter H{ oo ) is 
given, we can retrieve the pulse amplitude and phase very rapidly without using 
iterative calculations. Such experiment was actually demonstrated and will be 
discussed later. 

Ill LIMIT OF SONOGRAM CHARACTERIZATION OF OPTICAL PULSES 
[A] REQUIREMENT FOR THE FILTERING BANDWIDTH 
We discuss requirements for the bandwidth of the filter based on the pulse 
reconstruction formula (9). 
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We assume a chirped Gaussian pulse for the pulse under test: 

a(t) = exp[-£(l+jC)j , (10) 

where the time is normalized to the pulse width parameter, and C is the chirp parameter. 
We also assume the transfer function of the filter having a Gaussian distribution: 

"<•") = «p (") 
By using these Gaussian functions and (1), (2) and (3), the real sonogram 
G(oo,t) and its characteristic function M(6 , x) can be expressed in the following 
analytical forms: 

G(«,t) = exp rrrnJTrr X «P' 



i^ + iy + c* J (wg + i) J + c 2 J 

+ ly + c 2 J - 



(12) 



M(^) - ex P l- - 2 ^f + ^ |xexp[-^|xex P [-g^|. (13) 



Noting that 

^y = exp[-m±^]xe, P | p| r | 

we easily find that (9) gives the original pulse. 

It should be stressed that the co o-dependence of the sonogram and its characteristic 
function is cancelled out in (14), and the pulse waveform and phase, which are not 
dependent on coo, are retrieved. However, when co q < 1, the first term of (12) 
approaches to exp(- co o 2 t 2 ). This fact means that the temporal width of the sonogram is 
almost determined from the inverse of the filter bandwidth, and that the intrinsic 
information about the original pulse is masked by it. In such a case, the accurate pulse 
reconstruction becomes difficult since in the third term of (13), the factor of exp(- Q 
2 /4), which is necessary for the pulse reconstruction, is much smaller than the factor of 
exp(- 8 2 /4 co o 2 ), which must be cancelled out in (14). On the other hand, when co o > 1, 
ICI, we can not obtain sufficiently high spectral resolution to characterize the sonogram. 
Note that in the first term of (13), the factor of exp[- x 2 (l+C 2 )/4], which is essential to 
the pulse reconstruction, is much smaller than the factor of exp(- x 2 co o 2 /4), which must 
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be cancelled out in (14); hence, we can no longer retrieve the pulse accurately. We, 
thus, find the optimum value of co 0 — 1- 

[B] REQUIREMENT FOR THE SAMPLING PULSEWIDTH 
We next consider the requirement for the sampling pulse width. Let the sampling 
pulse have the Gaussian intensity waveform given by 

/.(*) = exp [-^] , (15) 

where 7s denotes the normalized pulse width parameter of the sampling pulse. 
The sonogram measured in Fig.l is given as 

G m (w,*) = J G{u,T)l,(T-t)dT. (16) 

The characteristic function Mm{ 6 , x ) of the measured sonogram Gm{ co , f) is given 
from (16) as 

M m .(9 1 T) = M(9.T)X a {9)\ (17) 
where T s( 9 ) denotes the Fourier transform of 7s(t). 
Using the Fourier transform of 7s(t) expressed as 

= «p \~\ - (18) 
the characteristic function for the measured sonogram is given from (3) and (17) as 

M,^ ) = e, P [- ^'+ : o' + O j xoxp [_«r] x ^[_fii±^5Ml]. (19) 

Then, we have 

M s (&,t) \ t 2 {l + C z )] f COt] f (i+T>)0 2 



A h( -e,t) = exp [ — J x cxp ["—J x exp [ • (2 °) 

Comparing (14) and (20), the requirement for reconstructing the pulse precisely is that 
Ts < 1 . This means that the sampling pulse width must be much shorter than the width 
of the pulse under test. However, when we know the sampling pulse shape and its 
Fourier transform in advance, we can deconvolute the measured characteristic function 
by using (17). This decon volution process is effective so long as Ts is comparable with 
or smaller than the pulse width under test. 

One may expect that when the bandwidth of the filter becomes narrower, the 
sonogram can be measured more precisely because the width of the filtered pulse 
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becomes wider than the sampling pulse width, allowing the pulse to be reconstructed. 
This statement is partially correct since the third term of (19) approaches to exp(- 6 2 /4 
co o 2 ) ,which is independent of 7s, as to 0 tends to zero. However, as mentioned before, 
this term does not contain the intrinsic information about the pulse, and is cancelled out 
in the pulse reconstruction process as shown in (20). Therefore, it has no meaning to use 
a filter bandwidth too small for the sonogram measurement and the succeeding pulse 
reconstruction. 



Reid deals with pulse retrieval from the sonogram measured in the 
cross-correlation setup shown in Fig.2, in which the sampling pulse is identical to the 
pulse under test. When the filter bandwidth is smaller than the spectral width of the 
pulse under test, the width of the frequency-filtered pulse usually becomes wider than 
the width of the pulse under test. Hence, it seems reasonable to expect that we can 
obtain the sonogram sufficiently accurate for pulse reconstruction. However, by 
following the method described in the previous subsection, we can show that pulse 
retrieval is not necessarily possible in this case. 

The sonogram measured in Fig 2 is given as 



where 7(t)=ls(r)l 2 is the intensity waveform of the pulse under test. Now, our problem is 
as follows: Can we really retrieve the pulse under test from Gm( co , t), instead of using 
G(co, 0? 

The characteristic function Mm( 6 , z) of the measured sonogram Gm( co , t) is 
given from (21) as 



where T ( 0 ) denotes the Fourier transform of I(t). For the Gaussian waveform given 
by (10), we have 



[C] SONOGRAM MEASUREMENT USING THE PULSE UNDER TEST 



AS THE SAMPLING PULSE 




(21) 



M m {e.r) = M(0.r)l(8)", 



J(0)=exp ^--J. 
Substitution of (13) and (23) into (22) yields 




(23) 
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M m (0 ; r) - exp ^ ^ J x exp ^ — j x exp I — . (24) 

We reconstruct the pulse from the measured sonogram Gm( co , t). Noting that 
M m {9,t) _ [ t 2 (l + C" 2 )] [ Cft] r B 2 



r cw] 

' 5- x exj 



" eXp [ ~1 1 xex P|-— I x ex P I-" I - - 

and substituting (25) into (9), we can obtain the complex amplitude of the reconstructed 
pulse as 

*>-«p[-^(i +J -c^)]. 

This reconstructed pulse is different from the pulse under test. Even if we use the 
iterative algorithm for pulse reconstruction assuming that the frequency window 
function is unknown, the retrieved pulse should be given by (26), which differs from 
the pulse under test. This result also denies the statement that we can measure the 
sonogram sufficiently for pulse reconstruction when the filter bandwidth is smaller than 
the spectral width of the pulse under test. We may apparently obtain an accurate 
sonogram by narrowing the filter bandwidth, but the original pulse waveform and phase 
are not reconstructed as already explained in the previous subsection. 

However, we can obtain the reconstructed pulse closer to the pulse under test, 
modifying the reconstruction process as follows. Fig. 3 shows the block diagram of 
such a process. As a first step, we use M 0 ( 8 , r ) = Mm( 6 , x ) for pulse reconstruction 
using (9). We next calculate 7 0 (t) and T Q (d) from the reconstructed pulse s Q (t) . The 
characteristic function is then modified as 

Using the modified characteristic function Mi and (9), we obtain the pulse si(t). As 
shown in Fig. 3, this process is repeated until the converged pulse is obtained. 

We apply this modification process to the chirped Gaussian pulse. When ICI CI, 
this process is very effective, and Table I shows the pulse width and chirp parameters of 
the reconstructed Gaussian pulse, which are normalized to the original values, as a 
function of the number of iteration. We find that these parameters rapidly converge at 
the real values. Even when we do not apply the process (the number of iteration^), the 
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ratio of the reconstructed pulse width to the original value is J~4/3, and the error is as 
small as 15 %. 

However, when 1CI > 1 , the reconstructed pulse moves toward 

s(t) = exp{-— -J—J. (28) 

and the actual pulse is no longer reconstructed. 

We, thus, conclude that only the pulse whose chirp parameter is small enough can 
be reconstructed from the sonogram measured in Fig.2. 

Table I 



number of iteration 


pulse-width parameter 


chirp parameter 
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IV. OPTICAL SAMPLING SYSTEM HAVING THE FUNCTION 
OF SONOGRAM CHARACTERIZATION 
In future ultra-high-speed optical fiber communication systems employing optical 
time-division multiplexing (OTDM), picosecond or sub-picosecond optical pulses will 
be transmitted. In these systems, the dispersive effect of optical devices such as fibers 
for transmission and optical filters induces serious waveform distortion and chirp of 
transmitted pulses. 

In order to diagnose the intensity waveform of such optical pulses, the optical 
sampling system is the most powerful tool, provided that we can prepare a sampling 
pulse, which has a width narrower than that of the pulse under test and is 
highly-synchronized with the pulse under test. 

On the other hand, there is strong demand for chirp measurement of optical pulses, 
and one of the methods to meet this demand is the sonogram characterization of optical 



11 



pulses. In all of the previous reports, after a pulse under test is frequency-filtered, the 
intensity waveform of the filtered pulse, which is called the sonogram, is measured by 
cross-correlating the filtered pulse with the original pulse under test. However, the 
actual sonogram cannot be obtained by this method, because the temporal resolution is 
limited by the shape of the pulse under test. One the contrary, when a short sampling 
pulse synchronized with the pulse under test is available, as is the case of the optical 
sampling system, we can realize precise sonogram characterization of the pulse under 
test by using the sampling pulse. 

Provided that the sonogram is measured by using experimental setup shown in 
Fig.4A, and that the sampling pulse width is much shorter than that of the pulse under 
test, we can measure the sonogram precisely. The pulse under test is completely 
retrieved from the measured sonogram by using (9). 

This experimental setup is regarded as an optical sampling system including the 
function of sonogram characterization, and can easily produce the following modified 
versions. In Fig. 4A, the sampling pulse is obtained by compressing the pulse under test 
itself. In Fig.4B, the ultrashort sampling pulse is incident on an optical device under 
test (DUT). Since the sonogram of the broadened output pulse is measured by 
cross-correlation using the sampling pulse, we can determine the impulse response of 
the DUT. 

V. OPTICAL SAMPLING SYSTEM AT 1.55. u m FOR THE MEASUREMENT 
OF PULSE WAVEFORM AND PHASE 
EMPLOYING SONOGRAM CHARACTERIZATION 
Practical implementation of such an optical sampling system at 1.55 fi m having 
the sonogram characterization function will be described. We demonstrate the 
measurement of impulse response of an optical bandpass filter as a specific application 
of the system. In the experimental setup, we first prepare a 200-fs optical pulse. Such 
pulse is incident on an optical bandpass filter under test, and the sonogram of the output 
pulse is measured by a highly-sensitive optical cross-correlator using two-photon 
absorption (TPA) in a Si avalanche photodiode (APD). The 200-fs pulse is used as a 
sampling pulse in the cross-correlator. The intensity and phase of the output pulse are 
very rapidly reconstructed from the sonogram by using a newly derived pulse 
reconstruction formula, enabling us to characterize the impulse response of the filter. 
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[A] SONOGRAM MEASUREMENT 



Fig. 5 shows the experimental setup. A Fourier-transform-limited 200-fs pulse 
having a 10-GHz repetition rate, a center wavelength of 1550 run, and a spectral width 
of 20.6 nm was obtained by supercontinuum compression of a mode-locked 
semiconductor-laser pulse. This pulse was branched into two paths. In one of the paths, 
an adjustable time delay was inserted, and the output from this paths was used as a 
sampling pulse was -10 dBm. In the other path, a three-cavity optical bandpass filter 
under test was inserted. The center wavelength of the filter was 1554.5 nm and the 3-dB 
bandwidth was 1.25 nm. The output pulse was amplified up to the average power of 10 
dBm, and incident on a tunable bandpass filter (BPF) with a 1-nm bandwidth for 
frequency gating. Fig. 6 shows the measured intensity and phase responses of the filter. 
The sampling pulse and the frequency-filtered pulse were combined and led to a 
cross-correlator using two-photon absorption in a Si APD. 

The sonogram trace was measured by sweeping the center frequency of the 
frequency gate and the delay time. The data points taken in such measurement were 256 
X 256. Fig. 7 shows the measured sonogram trace. Solid curves represent contours, 
where the normalized intensity is 0.2, 0.4, 0.6, 0.8 and 1. 



We first derive a pulse reconstruction formula from the sonogram closely following 
the method given by L.Cohen, "Time -frequency distributions- A review," Proc. IEEE, 
vol.77, no.7, pp.941-981, 1989. Let the complex amplitude of the signal pulse under 
test be S( 9 ) in the frequency domain and the complex transfer function of the filter be 
H( 0 ). When the center frequency of the filter is co , the sonogram P(t, co ) is given as 



The signal pulse under test in the frequency domain can be obtained from the following 
formula: 



[B] PULSE RECONSTRUCTION PROCESS 




(29) 




(30) 



where the M( 0 , z ) and A H ( 6 , t ) are defined as 




(31) 
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and 

A H (6, r) = J H* [u + 0 H (u - 0 e*"du. (32) 

If the transfer function of the filter A H ( 6 , r ) is known, the signal pulse under test 
can be reconstructed from the measured sonogram P(t, co ) using Eq.(30) without 
iterative calculations. The signal pulse in the time domain can be obtained with the 
inverse Fourier transformation of S( 8 ). Note that time-consuming iterative 
calculations have been indispensable for pulse reconstruction from the sonogram in the 
algorithm proposed in Reid. 

The pulse output from the filter under test was reconstructed from the measured 
sonogram (Fig.7) and the transfer function of the filter (Fig.6) by using Eq.(30). Figure 
4 shows the intensity waveform and phase of the reconstructed pulse. Since iterative 
calculations are not necessary, computation time for pulse reconstruction was shorter 
than 1 s, assuming a base 800 MHz Pentium ®m. 

The intensity waveform has an oscillatory structure in the leading edge. On the 
other hand, the phase response in the leading edge has abrupt tz -rad shifts when the 
intensity becomes zero. These characteristics clearly show the effect of the negative 
dispersion slope ( j3 3 < 0) of the bandpass filter. 

The auto-correlation trace calculated from the reconstructed pulse is shown in 
Fig.9A by circles, whereas the solid curve in Fig.9B is the directly measured one. 
Agreement between them is very good. On the other hand, the spectrum calculated 
from the reconstructed pulse is shown in Fig. 9B by circles. The solid curve represents 
the directly measured spectrum, where the fine structures corresponds to the 10-GHz 
repetition rate of the pulse train. The spectrum of the single pulse is given by its 
envelope, which is in good agreement with the circles. From these results, we find that 
the signal pulse is precisely reconstructed from its sonogram. 

We have constructed an optical sampling system at 1 .55 jj.m which enables us to 
measure the pulse waveform and phase through sonogram characterization. The 
measurement of impulse response of an optical bandpass filter is actually demonstrated 
by using this system. In our system, a 200-fs optical pulse is incident on an optical 
bandpass filter under test. The sonogram of the output pulse is measured by an optical 
cross-correlator using two-photon absorption in a Si avalanche photodiode, in which 
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the 200-fs pulse is also used as a sampling pulse. The intensity and phase of the output 
pulse are very rapidly reconstructed from the sonogram by using a newly derived pulse 
reconstruction formula. The measured intensity and phase responses clearly show the 
effect of the negative dispersion slope of the filter. 

While the preferred embodiment of the invention has been illustrated and 
described, it will be appreciated that various changes can be made therein without 
departing from the spirit and scope of the invention. 



